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Abstract 


The mitochondrial cytochrome c oxidase subunit 1 (COI) nucleotide sequences of Unaspis mabilis Lit & Barbecho 
(Hemiptera: Diaspididae), are provided for the first time. The total genomic DNA of each scale insect was extracted from 
individuals infesting lanzones leaves from three selected sites in Los Baños, Laguna. A partial COI gene amplicon with 
approximately 750 bp was obtained using the primer pair PcoF1 and LepR1. Nucleotide sequence alignment showed no 
variation among the COI sequences from all the samples. BLASTn search yielded no significant match with any of the 
available sequences for Unaspis species. The closest hit was Aulacaspis tubercularis Newstead (GenBank Accession No. 
HM474091) with 87.4% nucleotide similarity. Nonetheless, phylogenetic analyses revealed that generated COI sequences 
from U. mabilis form a monophyletic clade with U. yanonensis and U. euonymi, with closer proximity to the former. 
These findings also strengthen the species status of U. mabilis under the genus Unaspis. The DNA barcodes generated 
from this study (GenBank Acc. Nos. MN114099, MN14101, and MN114102), could, therefore, be used to verify the 
species identity of other lanzones scale accessions, as well as monitor the distribution and spread of U. mabilis, which 


would greatly influence possible pest management options. 
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Introduction 


Lanzones, Lansium parasiticum (Osbeck) Sahni & Bennet 
(=Lansium domesticum Corréa), is a well-known fruit sold in 
the Philippine market. In 2013, the volume of production for 
lanzones reached approximately 35,000 MT (Philippine 
Statistics Authority 2018). Aside from being a table fruit, it is 
also known to have medicinal value, which makes it a good 
choice of crop for massive production (Tilaar et al. 2008). 
Lanzones is known as an alternative medicine to treat illnesses 
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such as diarrhea, fever, dysentery, and even malaria (Yapp and 
Yap 2003). However, recent threats posed by insect pests 
gravely affected its local production. 

The recently described Unaspis mabilis Lit & Barbecho 
(Hemiptera: Diaspididae) is a species of armored scale insect 
that infests on lanzones. Serious infestation was first observed in 
2004 in North Cotabato (Provido 2007). The same scale insect 
was also found to attack lanzones in Davao (Provido 2007; DA- 
SMIARC 2014), South Cotabato and Saranggani provinces 
(Tacio 2008). In addition, an outbreak in Makilala, North 
Cotabato was reported in 2009 (Lit & Barbecho 2014). The 
infestation was reported to have spread to Negros Occidental 
(pers. comm.), Batangas (Watson 2015), Laguna (Lit & 
Barbecho 2014; Watson 2015), and recently in Aklan (Lit 2015, 
unpublished identification report). 

Owing to the minute nature of these scale insects, there 
were several misidentifications, as summarized by Lit & 
Barbecho (2014) and Watson (2015). They were initially 
misidentified in popular pest reports as Lepidosaphes ulmi 
(Linnaeus), and later as Lepidosaphes beckii (Newman), as their 
scale cover superficially resembles those of mussel scales (e.g. 
Mindanews 2007; Provido 2007; Cahatian & Pagal 2011). Later 
on, it was thought to represent Unaspis citri (Comstock) (e.g. 
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BPI-NPQD 2014). There were also a few misconceptions in 
national news that “coconut scale insects Aspidiotus spp. 
spreads to lanzones” (Valencia 2014). However, a careful and 
thorough morphological examination by Lit & Barbecho (2014) 
established that the mussel scale insect is a new species, 
Unaspis mabilis Lit & Barbecho. 

Despite only being reported from the country, Lit & 
Barbecho (2014) hypothesized that it may be an invasive 
species, based on its aggressive behavior and spread. With that, 
a more robust monitoring system should be implemented. The 
use of a DNA barcode for U. mabilis, would serve as a useful 
tool to survey and monitor its spread as well as understand its 
biogeographical origins, which is essential in designing efficient 
management programs. 

The mitochondrial gene, cytochrome c oxidase I (COD, 
serves as a standard DNA barcode for animals. This COI gene 
has a simple genomic structure, is distinct, easy to isolate and 
has a rapid pace of evolutionary change (Avise et al. 1987), and 
is one of the most frequently sequenced mitochondrial DNA 
genes (Caterino et al. 2000). Here, the partial sequences of the 
mitochondrial COI gene that could be used as a DNA barcode 
of the lanzones scale insect, U. mabilis, is provided. 


Materials and Methods 


Insect samples 

Unaspis mabilis samples were collected from lanzones 
trees in three sites in Los Bafios, Laguna: Mt. Makiling Station 
9, Bagong Silang, and Mudspring. High rates of lanzones scale 
insect infestation were recorded in these areas. The specimens 
were verified as U. mabilis by Dr. Ireneo L. Lit, Jr., the main 
author of the species. 


Total genomic DNA extraction 

The total genomic DNA of each lanzones scale insect was 
extracted following the protocol of the Animal and Fungi 
Genomic DNA Extraction Kit™ (Jena Bioscience, Germany) 
with slight modifications. Briefly, fresh specimens were placed 
individually on the cap of a 1.5 mL microcentrifuge tube and 
homogenized in 40 uL of cell lysis buffer with 0.3 uL 
proteinase K solution. The homogenates were incubated at 55°C 
for 2 h. After which, 30 uL of protein precipitation solution was 
then added into each tube, mixed on a vortex mixer for 30 s and 
centrifuged at 15,000 x g for 10 min. The supernatant was then 
transferred into a new microcentrifuge tube with 80 uL of 
isopropanol, which was then mixed by inverting the tubes 50 
times. The genomic DNA was collected by centrifugation at 
15,000 x g for 15 min to form the DNA pellet. After discarding 
the supernatant, 100 ul washing buffer was added followed by 
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centrifugation at 15,000 x g for 20 min. The DNA pellet was 
allowed to air-dry and was resuspended in 20 uL DNA 
hydration solution with 0.3 uL of RNAse. The tubes were 
incubated at 37°C for 1 h and then incubated at 65°C for another 
hour. The samples were stored at -20°C until use. 


Amplification of COI gene 

The COI gene was amplified using PcoFl (5’- 
CCTTCAACTAATCATAAAAATATYAG-3’), and LepR1 (5’- 
TAAACTTCTGGATGTCCAAAAAATCA -3°) primers 
following the thermocycle conditions; 2 min at 95°C; 5 cycles of 
40 s at 94°C, 40 s at 45°C, 70 s at 72°C; 40 cycles of 40 s at 94° 
C, 40 s at 51°C, 70 s at 72°C; 5 min at 72°C; held at 4°C (Park 
et al. 2010) in 25 uL PCR reaction mixture containing 12.5 uL 
of 2X Taq Master mix (Vivantis, Malaysia), 10 pmol PcoF1, 10 
pmol LepR1, 4 mM MgCh, and 50 ng genomic DNA. The 
amplicons were resolved in 1% agarose gel in 0.5x TBE buffer 
stained with GelRed® (Biotium, Inc., Fremont, CA, USA). 
Molecular size of the PCR products were estimated using 
VC100 bp Plus DNA Ladder (Vivantis, Malaysia). The gels 
were viewed and photographed using Alpha Imager Mini 
(Protein Simple, San Jose, CA, USA) at 302 nm. 


COI nucleotide sequence analysis 

Samples with corresponding expected molecular weight of 
amplicons were sent to Apical Scientific Sdn. Bhd. (Taman 
Serdang Perdana, 43300 Seri Kembangan, Selangor, Malaysia) 
via Asia Gel Corp. (Villa Lourdes Townhomes, Congressional 
Avenue, Barangay Bahay Toro, Quezon City 1100) for Sanger 
sequencing. The generated nucleotide sequences were edited in 
BioEdit software (Hall 1999). Those with high quality 
nucleotide sequences were then subjected to BLASTn (https:// 
blast.ncbi.nlm.nih.gov/Blast.cgi7PAGE _TYPE=Blast Search) 
for NCBI database mining. Phylogenetic analysis of the 
generated COI sequences was also performed using Mega X 
(Kumar et al. 2018) 


Results and Discussion 


The primer pair PcoF1 and LepR1, designed to address the 
low recovery of COI barcodes for scale insects (Park et al. 
2010), successfully generated an approximately 750 bp COI 
gene amplicon from Unaspis mabilis samples from the three 
collection sites (Figure 1). The amplicons generated were of the 
same expected molecular size as those reported by Park et al. 
(2010). 

High-quality COZ nucleotide sequences of 649 bp in length 
were obtained from the insect samples from Barangay Bagong 
Silang (n=10), Mudspring Area (n=12), and Mt. Makiling 
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Figure 1. Molecular size estimation of the amplicons of cytochrome c 
oxidase I (COL) gene region of Unaspis mabilis Lit & Barbecho from 
three sites in Los Bafios, Laguna (Lanes 1-12 in each panel). The PCR 
products are resolved in 1% agarose gel ran in 0.5x TBE buffer and 
stained with GelRed® (Biotium, Inc., USA). VC 100 bp Plus DNA 


Ladder (Vivantis, Malaysia) was used as a molecular size standard 
(MW). 





Station 9 (n=10). Alignment of the nucleotide sequences 
showed a 100% identity among all the samples without any stop 
codon within the deduced amino acid sequence (Figure 2) and 
thus, represents an uninterrupted open reading frame. This 
result confirmed the integrity of the COI nucleotide sequences 
generated, and hence, could serve as an ideal DNA barcode for 
U. mabilis. 

Comparison of the U. mabilis COI nucleotide sequence 
with other COI sequences in GenBank revealed an 87.40% 
identity match with another armored scale, Aulacaspis 
tubercularis Newstead (GenBank Acc. No. HM474091.1). The 


relatively low nucleotide identity match with A. tubercularis is 
due to 86 substitution sites observed within the COJI gene region. 
Surprisingly, an even lower percentage nucleotide identity 
match (86.35%) was observed between the U. mabilis COI 
nucleotide sequences and its congener, U. euonymi (Comstock) 
(GenBank Acc. No. HM474405.1), with 91 substitution sites 
observed (Figure 3). These results further proved the utility of 
the generated COI sequences as DNA barcode for U. mabilis. 
Despite a deep sequence genetic divergence (>2%) and 
apparent BLAST mismatch, phylogenetic analysis of 
Diaspididae COZI nucleotide sequences revealed that the 
generated U. mabilis DNA barcodes form a monophyletic clade 
with U. yanonensis (Kuwana) and U. euonymi indicative of its 
valid species status under Unaspis genus (Figure 4). 
Interestingly, U. mabilis had a closer proximity to U. yanonensis 
than U. euonymi. Nevertheless, the degree of divergence and 
uniqueness of U. mabilis COI barcode from the rest of the 
sequences examined proved its utility as a marker for this 
species. In a DNA barcoding study using the same COI locus 
done by Park et al. (2011b), high genetic differences were also 
observed among scale insects belonging to Pseudococcidae and 
Diaspididae. It was reported that genetic divergences within 
species (n=57 taxa) and in congeneric species (n=30 taxa) 
ranged from 0-5.98% with a mean of 0.97%, and 1.93-23.29% 
with a mean of 10.09%, respectively. Almost similar degree of 
divergence ranges within and in congeneric species was also 
observed in 1,090 Heteroptera COI barcodes (Park et al. 201 1a). 
The nucleotide sequence identity matrix of U. mabilis, U. 
euonymi, and A. tubercularis showing the divergence among the 
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Figure 2. Nucleotide sequence alignment and the deduced amino acid sequence of cytochrome c oxidase I 
gene of Unaspis mabilis Lit & Barbecho infesting lanzones Lansium parasiticum collected in Bagong Silang 
(MN114099), Mud Spring (MN114101) and Mt. Makiling Station 9 (MN114102), Los Baños, Laguna. 
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Figure 3. Multiple nucleotide sequence alignment of cytochrome c oxidase I gene consensus 


sequence of Unaspis mabilis Lit & Barbecho infesting lanzones, Lansium parasiticum, Aulascaspis 
tubercularis (HM474091.1) and U. euonymi (HM474405.1). 


species is summarized in Table 1. U. mabilis samples are 
significantly different from U. euonymi and A. tubercularis with 
genetic distances of 14.4% and 13.3%, respectively. These 
sequence variation values are highly significant since a 
minimum of 2.0% divergence in the COI nucleotide sequence is 
enough for species delineation (Hebert et al. 2003). 

By itself, the closeness of U. mabilis with U. euonymi is 
attributable to their being congeners but the significant 
difference confirms the distinctness of U. mabilis as a species, 
as first established by Lit & Barbecho (2014) based on 
morphological characters. In fact, it is the first Unaspis species 
to be reported on a host plant in the Meliaceae. The closest 
match, albeit 87% only, with A. tubercularis, can be attributed 


to the closeness of the subtribe Chionaspidina (where 


Aulacaspis belongs) and the Unaspis group, within the tribe 
Diaspidini (Normark et al. 2019). 

With the very minute size of U. mabilis, the availability of 
a DNA barcode will speed up studies on its distribution across 
geographical locations. This DNA barcode, which is the first 
report for U. mabilis, could also be used as a tool to monitor and 
trace the spread of this insect as well as identify other potential 
host plants. The identical COZ sequences obtained from the three 
populations suggest the presence of only one haplotype, an 
indication that the establishment is recent. Comparison of the 
DNA barcode from U. mabilis from other locations could, 
therefore, provide potential data for understanding its origins. 
Such information is essential in designing pest management 
strategies against this pest. 


Table 1. Nucleotide sequence identity matrix of Unaspis mabilis Lit & Barbecho infesting lanzones, Lansium parasiticum collected in 
Bagong Silang (MN114099), Mud Spring (MN114101) and Mt. Makiling Station 9 (MN114102), Los Baños, Laguna. 


U. mabilis 


COL Sequences MN114099  MN114101 
U. mabilis MN114099 1 
U. mabilis MN114101 1 
U. mabilis MN114102 1 1 
A. tubercularis HM74091.1 0.867 0.867 
U. euonymi HM74405.1 0.856 0.856 
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U. mabilis 


U. mabilis A. tubercularis U. euonymi 


MN114102 HM74091.1 HM74405.1 
1 0.867 0.856 
1 0.867 0.856 
0.867 0.856 
0.867 0.855 
0.856 0.855 
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Chionaspis pinifoliae KR043790.1 
Chionaspis pinifoliae KR042634.1 
Chionaspis pinifoliae KR041932.1 


Chionaspis pinifoliae KR042165.1 


Pinnaspis strachani HM474274.1 
Pinnaspis strachani GU936948.1 
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Lepidosaphes pini HM474218.1 
Lepidosaphes ulmi KX091211.1 
Unaspis euonymi HM474405.1 
Unaspis yanonensis HM474409.1 
Unaspis mabilis MN114099 
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Pseudaulacaspis nr. pentagona KX091218.1 
Pseudaulacaspis nr. pentagona KX091216.1 


Pseudococcus jackbeardsleyi KY373149.1 
Planococcus citri EU250560.1 


Figure 4. Phylogenetic relationship of Unaspis mabilis (n=3) and selected Diaspididae species (n=27) based on partial 
nucleotide sequences of cytochrome C oxidase I sequences analyzed by neighbor-joining using K2P distance model and 1000 
bootstraps, with Pseudococcidae species namely, Pseudococcus jackbeardsleyi (KY373149.1) and Planococcus citri 


(EU250560.1) as outgroups. 


Conclusion and Recommendation 


A 649 bp long COI gene sequence was generated as a 
DNA barcode for U. mabilis. This important molecular 
information can facilitate studies on the distribution of this 
species across geographical locations, particularly where 
lanzones are grown. It could also provide potential data for 
understanding its biogeographical origins, especially since there 
is doubt as to the species current endemic status in the 
Philippines and the hypothesis that it might be an invasive alien 
species. In turn, such understanding can have potentially great 
impacts on possible pest management options. Future detection 
on other potential host plants can also be hastened using this 
molecular barcode. 
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